Genetic ancestry and diagnostic yield of exome sequencing in a diverse population.

It has been suggested that diagnostic yield (DY) from Exome Sequencing (ES) may be lower among patients with non-European ancestries than those with European ancestry. We examined the association of DY with estimated continental/subcontinental genetic ancestry in a racially/ethnically diverse pediatric and prenatal clinical cohort. Cases (N = 845) with suspected genetic disorders underwent ES for diagnosis. Continental/subcontinental genetic ancestry proportions were estimated from the ES data. We compared the distribution of genetic ancestries in positive, negative, and inconclusive cases by Kolmogorov-Smirnov tests and linear associations of ancestry with DY by Cochran-Armitage trend tests. We observed no reduction in overall DY associated with any genetic ancestry (African, Native American, East Asian, European, Middle Eastern, South Asian). However, we observed a relative increase in proportion of autosomal recessive homozygous inheritance versus other inheritance patterns associated with Middle Eastern and South Asian ancestry, due to consanguinity. In this empirical study of ES for undiagnosed pediatric and prenatal genetic conditions, genetic ancestry was not associated with the likelihood of a positive diagnosis, supporting the equitable use of ES in diagnosis of previously undiagnosed but potentially Mendelian disorders across all ancestral populations.

SLCO1B1 functional variants and statin-induced myopathy in people with recent genealogical ancestors from Africa: a population-based real-world study.

Background: Clinical pharmacogenetic implementation guidelines for statin therapy are derived from evidence of primarily Eurocentric study populations. Functional SLCO1B1 variants that are rare in these study populations have not been investigated as a determinant of statin myotoxicity and are thus missing from guideline inclusion. Objective: Determine the relationship between candidate functional SLCO1B1 variants and statin-induced myopathy in people with recent genealogical ancestors from Africa. Design: Population-based pharmacogenetic study using real-world evidence from electronic health record-linked biobanks. Setting: Various health care settings. Participants: Self-identified white and Black statin users with genome-wide genotyping data available. Measurements: Primarily, the odds of statin-induced myopathy + rhabdomyolysis. Secondarily, total bilirubin levels. Thirdly, cell-based functional assay results. Results: Meta-analyses results demonstrated an increased risk of statin-induced myopathy + rhabdomyolysis with c.481+1G>T (odds ratio [OR] = 3.27, 95% confidence interval [CI] 1.43-7.46, P =.005) and c.1463G>C (OR = 2.45, 95% CI 1.04-5.78, P =.04) for Black participants. For White participants, c.521T>C was also significantly associated with increased risk of statin-induced myopathy + rhabdomyolysis (OR = 1.41, 95% CI 1.20-1.67, P =5.4x10 -5 ). This effect size for c.521T>C was similar in the Black participants, but did not meet the level of statistical significance (OR = 1.47, 95% CI 0.58-3.73, P =0.41). Supporting evidence using total bilirubin as an endogenous biomarker of SLCO1B1 function as well as from cell-based functional studies corroborated these findings. Limitations: Data limited to severe statin myotoxicity events. Conclusion: Our findings implicate Afrocentric SLCO1B1 variants on preemptive pharmacogenetic testing panels, which could have an instant impact on reducing the risk of statin-associated myotoxicity in historically excluded groups. Primary Funding Source: National Institutes of Health, Office of the Director - All of Us (OD-AoURP).

Diagnostic Yield of Exome Sequencing in a Diverse Pediatric and Prenatal Population is not Associated with Genetic Ancestry.

Purpose: It has been hypothesized that diagnostic yield (DY) from Exome Sequencing (ES) may be lower among patients with non-European ancestries than those with European ancestry. We examined the association of DY with estimated continental genetic ancestry in a racially/ethnically diverse pediatric and prenatal clinical cohort. Methods: Cases (N=845) with suspected genetic disorders underwent ES for diagnosis. Continental genetic ancestry proportions were estimated from the ES data. We compared the distribution of genetic ancestries in positive, negative, and inconclusive cases by Kolmogorov Smirnov tests and linear associations of ancestry with DY by Cochran-Armitage trend tests. Results: We observed no reduction in overall DY associated with any continental genetic ancestry (Africa, America, East Asia, Europe, Middle East, South Asia). However, we observed a relative increase in proportion of autosomal recessive homozygous inheritance versus other inheritance patterns associated with Middle Eastern and South Asian ancestry, due to consanguinity. Conclusions: In this empirical study of ES for undiagnosed pediatric and prenatal genetic conditions, genetic ancestry was not associated with the likelihood of a positive diagnosis, supporting the ethical and equitable use of ES in diagnosis of previously undiagnosed but potentially Mendelian disorders across all ancestral populations.

Diagnostic yield of pediatric and prenatal exome sequencing in a diverse population.

The diagnostic yield of exome sequencing (ES) has primarily been evaluated in individuals of European ancestry, with less focus on underrepresented minority (URM) and underserved (US) patients. We evaluated the diagnostic yield of ES in a cohort of predominantly US and URM pediatric and prenatal patients suspected to have a genetic disorder. Eligible pediatric patients had multiple congenital anomalies and/or neurocognitive disabilities and prenatal patients had one or more structural anomalies, disorders of fetal growth, or fetal effusions. URM and US patients were prioritized for enrollment and underwent ES at a single academic center. We identified definitive positive or probable positive results in 201/845 (23.8%) patients, with a significantly higher diagnostic rate in pediatric (26.7%) compared to prenatal patients (19.0%) (P = 0.01). For both pediatric and prenatal patients, the diagnostic yield and frequency of inconclusive findings did not differ significantly between URM and non-URM patients or between patients with US status and those without US status. Our results demonstrate a similar diagnostic yield of ES between prenatal and pediatric URM/US patients and non-URM/US patients for positive and inconclusive results. These data support the use of ES to identify clinically relevant variants in patients from diverse populations.

Lessons learned and recommendations for data coordination in collaborative research: The CSER consortium experience.

Integrating data across heterogeneous research environments is a key challenge in multi-site, collaborative research projects. While it is important to allow for natural variation in data collection protocols across research sites, it is also important to achieve interoperability between datasets in order to reap the full benefits of collaborative work. However, there are few standards to guide the data coordination process from project conception to completion. In this paper, we describe the experiences of the Clinical Sequence Evidence-Generating Research (CSER) consortium Data Coordinating Center (DCC), which coordinated harmonized survey and genomic sequencing data from seven clinical research sites from 2020 to 2022. Using input from multiple consortium working groups and from CSER leadership, we first identify 14 lessons learned from CSER in the categories of communication, harmonization, informatics, compliance, and analytics. We then distill these lessons learned into 11 recommendations for future research consortia in the areas of planning, communication, informatics, and analytics. We recommend that planning and budgeting for data coordination activities occur as early as possible during consortium conceptualization and development to minimize downstream complications. We also find that clear, reciprocal, and continuous communication between consortium stakeholders and the DCC is equally important to maintaining a secure and centralized informatics ecosystem for pooling data. Finally, we discuss the importance of actively interrogating current approaches to data governance, particularly for research studies that straddle the research-clinical divide.

Predicting genes from phenotypes using human phenotype ontology (HPO) terms.

The interpretation of genomic variants following whole exome sequencing (WES) can be aided using human phenotype ontology (HPO) terms to standardize clinical features and predict causative genes. We performed WES on 453 patients diagnosed prior to 18 years of age and identified 114 pathogenic (P) or likely pathogenic (LP) variants in 112 patients. We utilized PhenoDB to extract HPO terms from provider notes and then used Phen2Gene to generate a gene score and gene ranking from each list of HPO terms. We assigned Phen2Gene gene rankings to 6 rank classes, with class 1 covering raw gene rankings of 1 to 10 and class 2 covering rankings from 11 to 50 out of a total of 17,126 possible gene rankings. Phen2Gene ranked causative genes into rank class 1 or 2 in 27.7% of cases and the genes in rank class 1 were all associated with well-characterized phenotypes. We found significant associations between the gene score and the number of years, since the gene was first published, the number of HPO terms with an hierarchical depth greater or equal to 11, and the number of Online Mendelian Inheritance in Man terms associated with the phenotype and gene. We conclude that genes associated with recognizable phenotypes and terms deep in the HPO hierarchy have the best chance of producing a high gene score and ranking in class 1 to 2 using Phen2Gene software with HPO terms. Clinicians and laboratory staff should consider these results when HPO terms are employed to prioritize candidate genes.

Investigating rare pathogenic/likely pathogenic exonic variation in bipolar disorder.

Bipolar disorder (BD) is a serious mental illness with substantial common variant heritability. However, the role of rare coding variation in BD is not well established. We examined the protein-coding (exonic) sequences of 3,987 unrelated individuals with BD and 5,322 controls of predominantly European ancestry across four cohorts from the Bipolar Sequencing Consortium (BSC). We assessed the burden of rare, protein-altering, single nucleotide variants classified as pathogenic or likely pathogenic (P-LP) both exome-wide and within several groups of genes with phenotypic or biologic plausibility in BD. While we observed an increased burden of rare coding P-LP variants within 165 genes identified as BD GWAS regions in 3,987 BD cases (meta-analysis OR = 1.9, 95% CI = 1.3-2.8, one-sided p = 6.0 × 10-4), this enrichment did not replicate in an additional 9,929 BD cases and 14,018 controls (OR = 0.9, one-side p = 0.70). Although BD shares common variant heritability with schizophrenia, in the BSC sample we did not observe a significant enrichment of P-LP variants in SCZ GWAS genes, in two classes of neuronal synaptic genes (RBFOX2 and FMRP) associated with SCZ or in loss-of-function intolerant genes. In this study, the largest analysis of exonic variation in BD, individuals with BD do not carry a replicable enrichment of rare P-LP variants across the exome or in any of several groups of genes with biologic plausibility. Moreover, despite a strong shared susceptibility between BD and SCZ through common genetic variation, we do not observe an association between BD risk and rare P-LP coding variants in genes known to modulate risk for SCZ.

A Large-Scale Association Study Detects Novel Rare Variants, Risk Genes, Functional Elements, and Polygenic Architecture of Prostate Cancer Susceptibility.

To identify rare variants associated with prostate cancer susceptibility and better characterize the mechanisms and cumulative disease risk associated with common risk variants, we conducted an integrated study of prostate cancer genetic etiology in two cohorts using custom genotyping microarrays, large imputation reference panels, and functional annotation approaches. Specifically, 11,984 men (6,196 prostate cancer cases and 5,788 controls) of European ancestry from Northern California Kaiser Permanente were genotyped and meta-analyzed with 196,269 men of European ancestry (7,917 prostate cancer cases and 188,352 controls) from the UK Biobank. Three novel loci, including two rare variants (European ancestry minor allele frequency < 0.01, at 3p21.31 and 8p12), were significant genome wide in a meta-analysis. Gene-based rare variant tests implicated a known prostate cancer gene (HOXB13), as well as a novel candidate gene (ILDR1), which encodes a receptor highly expressed in prostate tissue and is related to the B7/CD28 family of T-cell immune checkpoint markers. Haplotypic patterns of long-range linkage disequilibrium were observed for rare genetic variants at HOXB13 and other loci, reflecting their evolutionary history. In addition, a polygenic risk score (PRS) of 188 prostate cancer variants was strongly associated with risk (90th vs. 40th-60th percentile OR = 2.62, P = 2.55 × 10-191). Many of the 188 variants exhibited functional signatures of gene expression regulation or transcription factor binding, including a 6-fold difference in log-probability of androgen receptor binding at the variant rs2680708 (17q22). Rare variant and PRS associations, with concomitant functional interpretation of risk mechanisms, can help clarify the full genetic architecture of prostate cancer and other complex traits. SIGNIFICANCE: This study maps the biological relationships between diverse risk factors for prostate cancer, integrating different functional datasets to interpret and model genome-wide data from over 200,000 men with and without prostate cancer.See related commentary by Lachance, p. 1637.

Genetic ancestry, skin pigmentation, and the risk of cutaneous squamous cell carcinoma in Hispanic/Latino and non-Hispanic white populations.

Although cutaneous squamous cell carcinoma (cSCC) is one of the most common malignancies in individuals of European ancestry, the incidence of cSCC in Hispanic/Latinos is also increasing. cSCC has both a genetic and environmental etiology. Here, we examine the role of genetic ancestry, skin pigmentation, and sun exposure in Hispanic/Latinos and non-Hispanic whites on cSCC risk. We observe an increased cSCC risk with greater European ancestry (P = 1.27 × 10-42) within Hispanic/Latinos and with greater northern (P = 2.38 × 10-65) and western (P = 2.28 × 10-49) European ancestry within non-Hispanic whites. These associations are significantly, but not completely, attenuated after considering skin pigmentation-associated loci, history of actinic keratosis, and sun-protected versus sun-exposed anatomical sites. We also report an association of the well-known pigment variant Ala111Thr (rs1426654) at SLC24A5 with cSCC in Hispanic/Latinos. These findings demonstrate a strong correlation of northwestern European genetic ancestry with cSCC risk in both Hispanic/Latinos and non-Hispanic whites, largely but not entirely mediated through its impact on skin pigmentation.

The role of exome sequencing in newborn screening for inborn errors of metabolism.

Public health newborn screening (NBS) programs provide population-scale ascertainment of rare, treatable conditions that require urgent intervention. Tandem mass spectrometry (MS/MS) is currently used to screen newborns for a panel of rare inborn errors of metabolism (IEMs)1-4. The NBSeq project evaluated whole-exome sequencing (WES) as an innovative methodology for NBS. We obtained archived residual dried blood spots and data for nearly all IEM cases from the 4.5 million infants born in California between mid-2005 and 2013 and from some infants who screened positive by MS/MS, but were unaffected upon follow-up testing. WES had an overall sensitivity of 88% and specificity of 98.4%, compared to 99.0% and 99.8%, respectively for MS/MS, although effectiveness varied among individual IEMs. Thus, WES alone was insufficiently sensitive or specific to be a primary screen for most NBS IEMs. However, as a secondary test for infants with abnormal MS/MS screens, WES could reduce false-positive results, facilitate timely case resolution and in some instances even suggest more appropriate or specific diagnosis than that initially obtained. This study represents the largest, to date, sequencing effort of an entire population of IEM-affected cases, allowing unbiased assessment of current capabilities of WES as a tool for population screening.

A multiethnic genome-wide analysis of 44,039 individuals identifies 41 new loci associated with central corneal thickness.

Central corneal thickness (CCT) is one of the most heritable human traits, with broad-sense heritability estimates ranging between 0.68 to 0.95. Despite the high heritability and numerous previous association studies, only 8.5% of CCT variance is currently explained. Here, we report the results of a multiethnic meta-analysis of available genome-wide association studies in which we find association between CCT and 98 genomic loci, of which 41 are novel. Among these loci, 20 were significantly associated with keratoconus, and one (RAPSN rs3740685) was significantly associated with glaucoma after Bonferroni correction. Two-sample Mendelian randomization analysis suggests that thinner CCT does not causally increase the risk of primary open-angle glaucoma. This large CCT study explains up to 14.2% of CCT variance and increases substantially our understanding of the etiology of CCT variation. This may open new avenues of investigation into human ocular traits and their relationship to the risk of vision disorders.

Telomere length and socioeconomic status at neighborhood and individual levels among 80,000 adults in the Genetic Epidemiology Research on Adult Health and Aging cohort.

BACKGROUND: Telomere length (TL) may serve as a biologic marker of aging. We examined neighborhood and individual-level socioeconomic status (SES) in relation to TL. METHODS: The study included 84,996 non-Hispanic white subjects from the Genetic Epidemiology Research on Adult Health and Aging (GERA) cohort, part of the Research Program on Genes, Environment and Health. Relative TL (T/S) was log2 transformed to improve normality and standardized to have mean 0 and variance 1. Neighborhood SES was measured using the Neighborhood Deprivation Index (NDI), and individual SES was measured by self-reported education level. We fit linear regression models of TL on age, sex, smoking, body mass index, comorbidities, NDI, and education level. We tested for differences in the associations by sex and nonlinearity in the association of NDI with TL. RESULTS: Each SD increase in NDI was associated with a decrease of 0.0192 in standardized TL, 95% confidence interval (CI) = -0.0306, -0.0078. There was no evidence of nonlinearity in the association of NDI with TL. We further found that less than high school education was associated with a decrease of 0.1371 in standardized TL, 95% CI = -0.1919, -0.0823 as compared to a college education. There were no differences in the associations by sex. CONCLUSIONS: We found evidence that both lower neighborhood SES and lower individual-level SES are associated with shorter TL among non-Hispanic whites. Our findings suggest that socioeconomic factors may influence aging by contributing to shorter TL.

A Large Multiethnic Genome-Wide Association Study of Adult Body Mass Index Identifies Novel Loci.

Body mass index (BMI), a proxy measure for obesity, is determined by both environmental (including ethnicity, age, and sex) and genetic factors, with > 400 BMI-associated loci identified to date. However, the impact, interplay, and underlying biological mechanisms among BMI, environment, genetics, and ancestry are not completely understood. To further examine these relationships, we utilized 427,509 calendar year-averaged BMI measurements from 100,418 adults from the single large multiethnic Genetic Epidemiology Research on Adult Health and Aging (GERA) cohort. We observed substantial independent ancestry and nationality differences, including ancestry principal component interactions and nonlinear effects. To increase the list of BMI-associated variants before assessing other differences, we conducted a genome-wide association study (GWAS) in GERA, with replication in the Genetic Investigation of Anthropomorphic Traits (GIANT) consortium combined with the UK Biobank (UKB), followed by GWAS in GERA combined with GIANT, with replication in the UKB. We discovered 30 novel independent BMI loci (P < 5.0 × 10-8) that replicated. We then assessed the proportion of BMI variance explained by sex in the UKB using previously identified loci compared to previously and newly identified loci and found slight increases: from 3.0 to 3.3% for males and from 2.7 to 3.0% for females. Further, the variance explained by previously and newly identified variants decreased with increasing age in the GERA and UKB cohorts, echoed in the variance explained by the entire genome, which also showed gene-age interaction effects. Finally, we conducted a tissue expression QTL enrichment analysis, which revealed that GWAS BMI-associated variants were enriched in the cerebellum, consistent with prior work in humans and mice.

A multiethnic genome-wide association study of primary open-angle glaucoma identifies novel risk loci.

Primary open-angle glaucoma (POAG) is a leading cause of irreversible vision loss, yet much of the genetic risk remains unaccounted for, especially in African-Americans who have a higher risk for developing POAG. We conduct a multiethnic genome-wide association study (GWAS) of POAG in the GERA cohort, with replication in the UK Biobank (UKB), and vice versa, GWAS in UKB with replication in GERA. We identify 24 loci (P < 5.0 × 10-8), including 14 novel, of which 9 replicate (near FMNL2, PDE7B, TMTC2, IKZF2, CADM2, DGKG, ANKH, EXOC2, and LMX1B). Functional studies support intraocular pressure-related influences of FMNL2 and LMX1B, with certain Lmx1b mutations causing high IOP and glaucoma resembling POAG in mice. The newly identified loci increase the proportion of variance explained in each GERA race/ethnicity group, with the largest gain in African-Americans (0.5-3.1%). A meta-analysis combining GERA and UKB identifies 24 additional loci. Our study provides important insights into glaucoma pathogenesis.

A large electronic-health-record-based genome-wide study of serum lipids.

A genome-wide association study (GWAS) of 94,674 ancestrally diverse Kaiser Permanente members using 478,866 longitudinal electronic health record (EHR)-derived measurements for untreated serum lipid levels empowered multiple new findings: 121 new SNP associations (46 primary, 15 conditional, and 60 in meta-analysis with Global Lipids Genetic Consortium data); an increase of 33-42% in variance explained with multiple measurements; sex differences in genetic impact (greater impact in females for LDL, HDL, and total cholesterol and the opposite for triglycerides); differences in variance explained among non-Hispanic whites, Latinos, African Americans, and East Asians; genetic dominance and epistatic interaction, with strong evidence for both at the ABO and FUT2 genes for LDL; and tissue-specific enrichment of GWAS-associated SNPs among liver, adipose, and pancreas eQTLs. Using EHR pharmacy data, both LDL and triglyceride genetic risk scores (477 SNPs) were strongly predictive of age at initiation of lipid-lowering treatment. These findings highlight the value of longitudinal EHRs for identifying new genetic features of cholesterol and lipoprotein metabolism with implications for lipid treatment and risk of coronary heart disease.

Whole exome and whole genome sequencing with dried blood spot DNA without whole genome amplification.

Newborn screening (NBS) for rare conditions is performed in all 50 states in the USA. We have partnered with the California Department of Public Health Genetic Disease Laboratory to determine whether sufficient DNA can be extracted from archived dried blood spots (DBS) for next-generation sequencing in the hopes that next-generation sequencing can play a role in NBS. We optimized the DNA extraction and sequencing library preparation protocols for residual infant DBS archived over 20 years ago and successfully obtained acceptable whole exome and whole genome sequencing data. This sequencing study using DBS DNA without whole genome amplification prior to sequencing library preparation provides evidence that properly stored residual newborn DBS are a satisfactory source of DNA for genetic studies.

A large multi-ethnic genome-wide association study identifies novel genetic loci for intraocular pressure.

Elevated intraocular pressure (IOP) is a major risk factor for glaucoma, a leading cause of blindness. IOP heritability has been estimated to up to 67%, and to date only 11 IOP loci have been reported, accounting for 1.5% of IOP variability. Here, we conduct a genome-wide association study of IOP in 69,756 untreated individuals of European, Latino, Asian, and African ancestry. Multiple longitudinal IOP measurements were collected through electronic health records and, in total, 356,987 measurements were included. We identify 47 genome-wide significant IOP-associated loci (P < 5 × 10-8); of the 40 novel loci, 14 replicate at Bonferroni significance in an external genome-wide association study analysis of 37,930 individuals of European and Asian descent. We further examine their effect on the risk of glaucoma within our discovery sample. Using longitudinal IOP measurements from electronic health records improves our power to identify new variants, which together explain 3.7% of IOP variation.

Pierpont syndrome associated with the p.Tyr446Cys missense mutation in TBL1XR1.

We present a 7-year old male with severe delays, hypotonia and dysmorphic features who had striking, deep palmar and plantar creases and pillowing of the soft tissues of the palms and soles. His facial features included a high anterior hairline, small eyes with narrowed palpebral fissures, a bulbous nasal tip with a short columella, and a large mouth with a thin upper vermilion, and small chin. He had a submucous cleft palate, bilateral cryptorchidism and hydronephrosis. Cranial imaging demonstrated an Arnold Chiari malformation that was also present in his maternal uncle by report. Exome sequencing revealed a de novo heterozygous sequence variant, p.Tyr446Cys, in TBL1XR1 that has previously been reported in six patients with Pierpont syndrome. This sequence variant occurs in the carboxy-terminal, WD40 domain of the protein. As TBL1XR1 is a critical component of the NCoR/SMRT co-repressor complex and the WD40 repeats are hypothesized to interact with histone H2B and H4, the mutation may impact protein interactions necessary for stabilizing the complex with chromatin. De novo missense and frameshift mutations and deletions involving TBL1XR1 have been described in patients with intellectual disability and autism, but without any of the dysmorphic findings or malformations associated with Pierpont syndrome, implying a mutation-specific mechanism for the pathogenicity of p.Tyr446Cys. Our case is the first individual with this mutation to have a submucous cleft palate and hydronephrosis, although his severe delays, hypotonia, dysmorphic findings and emerging scoliosis appear consistent with previous reports. His distinctive facial and digital features are further evidence that p.Tyr446Cys results in a clinically recognizable, syndromic form of intellectual disability in contrast to other TBL1XR1 mutations.